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storage suggest stress-dependence of mutation rate, which 
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Abstract 

Background: Acceleration of adaptation dynamics by 
stress-induced hypermutation lias been found experimen- 
tally. Evolved evolvability is a prominent explanation. 
We investigate a more generally applicable explanation 
by a physical constraint. 

Methods and Results: A generic thermodynamical 
analysis of genetic information storage obviates physical 
constraints on the integrity of genetic information. The 
capability to employ metabolic resources is found as a 
major determinant of mutation probability in stored ge- 
netic information. Incorporation into a non-recombinant, 
asexual adaptation toy model predicts cases of markedly 
accelerated adaptation, driven by a transient increase of 
mutation rate. No change in the mutation rate as a ge- 
netic trait is required. The mutation rate of one and the 
same genotype varies dependent on stress level. 

Implications: Stress-dependent mutation rates 
are physically necessary and challenge a condition- 
independent genotype to mutation rate mapping. This 
holds implications for evolutionary theory and pathogen 
and cancer evolution. 
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Introduction 

Rationale: Conceptually, life can be understood as 
an informed, open organization of non-equilibrium 
thermodynamical processes, as the first chapters of 
most biochemistry textbooks will yield. This is a 
potent concept to understand the informational and 
energetic dealings of an organism, be it metabolism 
[30j . motility [19], the origin of life and early evo- 
lution [m ISl [m SH [13, evolution in general (4Q1, 
or pathogen evolution [Ml |36l |9] to name but a icw 
applications. 

We investigate the thermodynamical underpin- 
nings of storage of genetic information. In the light 
of adaptation dynamics, this analysis indicates that 
the adaptation of simple organisms can well be ac- 
celerated by a transient phase of stress-induced hy- 
permutation. Stress-induced hypermutation is com- 
monly understood as an organismic response that de- 
veloped as a consequence of a specific selection sce- 
nario. While this is a credible explanation, our analy- 
sis obviates a generally applicable physical constraint 
on genetic information storage as an explanation of 
stress-induced hypermutation. 

Background Genetic integrity and evo- 
lutionary adaptation: Mutations are a poten- 
tially lethal threat to hereditary information. Conse- 
quently, elaborate mechanisms to counteract muta- 
tions and their effects have evolved. However, evo- 
lution of these very mechanisms and life in general 
requires mutations as a source of genetic variation. 
Let this fundamental ambivalence guide our brief re- 
view of the evolutionary role of mutator genotypes 
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and stress- induced hypermutation. 

In both cases organisms have been observed to 
exhibit an elevated mutation rate when faced with 
an adaptation challenge. This suggests an organis- 
mic ability to respond to selection pressure with in- 
creased mutation rates. Classically mutation and se- 
lection are understood as two entirely independent el- 
ements of the evolutionary adaptation process. Thus, 
the experimental findings seem to have bearing on 
the foundations of our understanding of evolution- 
ary theory. This is the main similarity between mu- 
tator genotypes and stress-induced hypermutation. 
We will spend the rest of the introduction to clar- 
ify the fundamental differences of both concepts. We 
hope this prevents associations of stress-induced hy- 
permutation with mutator genotype studies, which 
we found to frequently occur. 

Background Selection of mutator geno- 
types: Elevated mutation rates can be found by 
comparison between species, strains and even cells 
of the same population. A vast body of literature 
has discussed, and in our opinion mostly explained 
this finding. For a general review see [42]. Mutator 
genotypes arise in response to specific selection sce- 
narios that favor higher mutation rates [71 [35l l33] . 
High mutation rate genotypes outplay low mutation 
rate genotypes in scenarios of frequently changing se- 
lection criteria [3j and intricate mechanisms for effi- 
cient and effective evolvability have evolved [15l [16] . 
Most importantly, elevated mutation rates of muta- 
tor genotypes are a selectable genetic trait. In this 
case the mutation rate is determined to a specific 
value from the genome of an organism, and this pre- 
cise value is the result of evolutionary adaptation to 
specific environmental conditions. 

Background Stress-induced hypermuta- 
tion: Cairn's Adaptive Mutation experiment [l2] 
and several following studies have shown that ge- 
netic mutation rates in organisms of the exact same 
genotype can increase transiently under stress. This 
leads to accelerated adaptation to the present stress, 
which is followed by reduction of mutation rates to 
pre-stress level. While it is known that a change 
in selection criteria can lead to increased variation 
by broader spread of the population frequencies into 
existent genotypes, these findings suggests that al- 



tered selection criteria transiently increase the rate 
of occurrence of new genotypes. This questions 
the fundamental notion of independence of mutation 
rate from selective pressure as shown in the classic 
Luria-Delbrck experiment [31]. Quite the opposite, 
the common explanation of the Cairn's experiment 
is that evolution itself has awarded organisms with 
mechanisms to master evolution more efficiently and 
effectively [16l [El [291 [22l [28l [231 [20l [10]. The ob- 
served mutation rate increase is not caused by a ge- 
netic difference, but rather a dynamic organismic re- 
sponse to externally applied stress. Genetic alter- 
ations in cellular functions in general can also impose 
a stress, which in turn induces a secondary hyper- 
mutation response. Note that this response is not 
directly mediated by a genetic alteration in genes di- 
rectly associated with mutation rate. 

Methods 

Thermodynamical analysis of mutation 
suppression 

General approach: The genome stores an organ- 
ism's hereditary information, and is therefore con- 
ceptually not different from any other information 
storage system: (1) The information storage system 
consists of a material medium that can take on dif- 
ferent configurations. A specific configurations corre- 
sponds to a specific information stored in it. In DNA 
this is a specific sequence of base pairs. (2) Only rare 
and very specific configuration serve to store informa- 
tion in a useful way. For example, two DNA strands 
of the same length that contain the same informa- 
tion must at each position have the exact same bases. 
(3) These specific configurations are low entropy non- 
equilibrium states. Randomizing influences drive the 
system towards equilibrium, which corresponds to 
full loss of information. Only free energy expendi- 
ture can maintain the storage system close to its orig- 
inal configuration/information. DNA is constantly 
exposed to mutagenic (randomizing) infiuences, and 
only metabolically expensive repair mechanisms can 
maintain the original configuration. A deviation from 
the original configuration corresponds to genetic mu- 
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tation. More generally speaking, the coupling of in- 
formation storage and the second law is an example 
of the generalized Carnot principle . 

Common sense might already suggest that both an 
increase of the mutagenic influence or an inability to 
provide or utilize metabolic resources for mechanisms 
of mutation suppression/repair will lead to stronger 
mutagenic alteration of genetic information. Still, in 
the following we will construct a formal thermody- 
namic toy model to outline the optimum mutation 
suppression capability given a certain level of muta- 
genic influence and a certain level of metabolic re- 
source (free energy) utilized for safe-guarding of the 
genetic information. The analysis is based on the 
definitions of a non-equilibrium thermodynamics de- 
scription of Master equations as presented in |25j . 

Thermodynamical analysis 

Let a binary sequence (6„) £ {0, 1} with n = 1, . . . ,N 
digits represent the binary genetic information of a 
generic organism. stands for a digit in its original 
value, while 1 represents a digit whose value has been 
changed by mutation. 

First, we will look at how the primary mutagenic 
influence and the metabolic resources utilized for the 
suppression of and protection against mutations set 
the probability p of a single digit to contain its origi- 
nal value, a state represented by the value 0. The ge- 
netic information storage system is a compartment at 
temperature T, which is in contact with a heat bath 
at temperature T^xt- Text represents the strength of 
the primary mutagenic influence. Omnipresent mu- 
tagenic influences like molecular thermal collisions, 
background radiation or effects from reaction with 
surrounding chemical species warrant for T^xt > at 
all times under all conditions. Heat flows into the in- 
formation storage system by diffusion and is exported 
by the use of metabolic resources: 

^=c{Text-T)-T-S'e, (1) 

where c is a coefficient characteristic for the diffu- 
sion of heat from the heat bath into the information 
storage compartment. S'^ — {dS/dt)c is the entropy 



production arising from the utilization of metabolic 
resource for the export of heat from the information 
storage system. For constant S'^ a globally stable 
steady state will be attained: 

dQ^O^T*^^. (2) 

Now as we know the steady state temperature T* 
of the information storage system, we still have to 
describe the dynamics between the original state 
and the mutated 1 state of the single binary digit. It 
seems reasonable to assume that the molecular struc- 
ture of the information carrier is not a priori biased 
for or against any of both possible states, so we assign 
the same forward and backward temperature depen- 
dent reaction rate kT* . fc is a generic transition rate 
between the and 1 state. The probability to be in 
state then evolves according to 

^ = kT*{l-2p), (3) 

which has a globally stable steady state at p* = 0.5. 
Let us further examine the thermodynamics of this 
system. The entropy has to be calculated for both 
states and 1: 

S{p) = \n{p) + ln(l -p) = \n{p - p^). (4) 
The time evolution of this entropy is 

S' = '4^'-l'i=kT*^l-^. (5) 
dt dp dt p — 

From the enumerator (1 — 2p)^ we see, that S" = 
at p* = 0.5. S" > everywhere else except p = and 
p — 1, where S' — > -|-oo. Thus, the entropy is a Lya- 
punov function on p, which shows a Lyapunov global 
stability of p* = 0.5 as well as = S{p*) = In 0.25. 
This corresponds to the second thermodynamical law 
of approach to equilibrium, or equivalently, maximum 
entropy of the system, when the system is not al- 
lowed to export its entropy. The S" — > +oo behavior 
at p — and p = 1 indicates that an infinite export 
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of entropy from the system would be required to sus- 
tain these states. These states also correspond to a 
perfectly determined value of the single binary digit, 
thus completely error-free storage is never possible. 

To maintain any useful information a probabil- 
ity p > p* — 0.5 is required. To provide this 
p > p* the storage system must be kept away from 
the p — p* equilibrium by a constant production of 
entropy S'^. (The index s stands for suppression of 
mutations.) The entropy export S'^ is achieved by 
utilization of metabolic resources. For constant S'^ 
a non-equilibrium steady state (NESS) is attained. 
The NESS p value can be calculated from insertion 
of the actual entropy export for mutation suppression 
S'^ into dSl): 



Substitution into ^ yields an explicit expres- 
sion for the optimal p given a certain utilization of 
metabolic resource S' 




C+(1 + S)2 



(9) 



where 's — S'^/c and C = IQkText/c. 

We insert the probability of a single binary digit 
not to be changed from its original value, i.e. not 
mutated, into an expression for the probability of m 
binary digits in a sequence of A'' to be mutated from 
their original value: 



where K — S'^/kT*. For ± we only consider the + 
alternative, which increases the probability of finding 
the original value of the digit. The — case would 
describe the case of an information storage system 
that utilizes metabolic resources to ensure mutation 
of digits. We will not bother to discuss this further 
in this work. 

Finally, let us assume that overall metabolic re- 
sources S'g for a single digit are available at a limited 
rate and can be allocated to cooling the information 
storage system by S'^ and suppressing mutations by 
Si: 

S'o — S'^ + S'g. (7) 

When we optimize this allocation for a maximal 
p value we can delineate the minimally possible per 
digit mutation rate for a given S'^. From dH) we see 
that a maximal p is attained when a maximal K value 
is taken. We use the above relation ([7]) together with 
the definition K = S'JkT* and ^ and optimize S'^ G 
(0,5;). We find the optimal S'^ = {S'^ - c)/2, which 
gives the maximal K: 

ckText ( 2 ) 



P^n-['l)[p{s')f~"'[i~p{s'T (10) 

This expression describes the probability of m mu- 
tations to occur in a sequence of N binary digits given 
our specific optimal assumptions. Fig. [T] shows an 
example for different Sq values. 

Stress-induced hypermutation in an 
adaptation toy model 

General approach 

In the last section we found a dependence of the mu- 
tation rate on the level of mutagenic influence and 
the availability and ability to make use of metabolic 
resources for mutation suppression. Now we need to 
understand what the consequences are on adaptation 
processes. We will assume no recombination, stasis, 
apoptosis etc. to keep the analysis straight forward. 
In reality these assumptions would be fulflUed by the 
simplest unicellular organisms. 

First, let us define how stress level and mutations 
in genetic information storage are connected: 

An elevated stress level a) reduces an organism's 
ability to deflect mutagenic influence and/or h) re- 
stricts or redistributes metabolic resources available 
to mutation suppression in an unstressed situation 
and/ or c) impairs the means to utilize these re- 
sources. Vice versa, any influence exerting one or 
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dS/dt 

Fig. 1: Probability of m mutations for dS/dt metabolic resource utilized on mutation suppres- 
sion: Probability Pm of m mutations in a binary sequence vs. entropy production S'^ ^ ^ for 
mutation suppression. The graph is an evaluation of (|10l) for sequence length iV = 200 with C = 200. 
The dashed black line represents no mutation m — 0. The heavy solid lines (Displayed in the follow- 
ing colors from right to left: gray, cyan, magenta, green, blue and red) represent the probabilities Pm 
of m = 1, 2, 3, 10, 20 and 30 mutations, respectively. 
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several of the aforementioned effects can he under- 
stood as stress. 

Adaptation toy model 

We are now equipped to construct a generic model of 
non-rccombinatory, mutational adaptation to a rapid 
change of living conditions. After constructing the 
general model we will execute and evaluate a numer- 
ical simulation. For a population of TV = 2^'' different 
binary genotypes g of sequence length N^. We ran- 
domly assign {dS/dt)g = S'g levels to each genotype 
g. S'g is the rate of entropy production from utiliza- 
tion of metabolic resource. We introduced before the 
assumption that S' is set by the stress an organism 
experiences under certain living conditions. Thus Sg 
can also be understood as a stress landscape. To in- 
vestigate two different adaptation scenarios we assign 
two types of stress landscape: m mild and r rough. 
We draw the values in two ways: 



(11) 



where Zg is drawn from a uniform distribiition 
P(0 < z < 1) = 1. This way both stress landscapes 
are distributed around the same mean value D. The 
mild one has a flatter distribution of S'g. The rough 
one has a few high peaks, and mostly low S'g. The 
mild case should model a gradual adaptation to a 
less harsh stress type. The rough case should model 
a sudden drastic change in living conditions with only 
few beneficial mutations having a strong adaptation 
effect. 

We avoid the commonly used fitness landscape be- 
cause it is a logical short circuit with respect to our 
introduced assumption: Fitness is defined as the rel- 
ative ability of an organism to have its genetic in- 
formation reproduced by the next generation. As we 
are trying to understand the effect of mutation rate 
differences on this very fitness, we should not define 
a reverse dependency of mutation rate on fitness. It 
seems conceptually easier to use the effective utiliza- 
tion of metabolic resource S' for mutation suppres- 
sion as a measure of the stress an organism is facing 
under a particular set of living conditions. 



Each genotype has a certain occupation number 
Ug (t) . This number can increase by infiow of all geno- 
types g G {1 . . ., N}, including g itself. We use a set 
of linear ordinary differential equations to model the 
temporal dynamics of the genotype occupation num- 
bers: 



^=Y^[GiS',)Pl^,^{S',)\. (12) 

k 

d{g, k) is the Hamming distance between binary 
genotypes g and k. The Hamming distance is the 
number of binary digits that are different between g 
and k. This is also the number of mutations by which 
g and k differ. G is the growth rate of genotype g, 
which we assume to be 



G{S') = 1-6 



-s' 



(13) 



We simulate and compare the simulation results for 
two types of populations: (1) A population with S'- 
dependent mutation rates, whose dynamics are de- 
scribed above (2) A "conventional" population, in 
which only the growth rates depend on S' , the mu- 
tation term is simply evaluated for the maximum 
max(S'g) of the whole stress landscape, no matter 
what the genotypes actual s'g is. 

We evaluate this toy model for A^g = 5 in a nu- 
merical simulation. The simulation is initialized with 
all ng{t = 0) = 0, only the one with the lowest S'g 
is assigned 1. This corresponds to a heavily destabi- 
lized wild type. Death and recombination effects are 
not considered. The dynamics were integrated with 
a simple first order Euler scheme, step size 0.1. 

The adaptation dynamics of population (1) and (2) 
are evaluated based on the mean S' =< S' >g across 
each population. Both populations are initiated and 
simulated with the same stress landscape S'g. We 
define a threshold for S', which is S' when the popu- 
lation (1) has increased half way from its initial t = 
value to its final value. For both populations we de- 
termine from the simiilation results the times T™2 
and (^^^ population (1) and (2), respectively), 

when this threshold is crossed. From these we can 
calculate a comparative value of the adaptation speed 
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-^1/2 

5Ti/2 > indicates faster adaptation of popula- 
tion (1), and STi/2 < indicates faster adaptation of 
population (2). 

Results 

Mutation rate depends on 
availability/utilization of metabolic 
resources 

As can be seen in Fig. [1] a decrease in the rate of en- 
tropy production by utilization of metabolic resources 
for protection of the hereditary information dS/dt 
leads to an increased probability to find higher num- 
bers of mutations in the genetic sequence. dS/ dt can 
be decreased for several reasons: The metabolic re- 
source is not available, the mechanisms to utilize the 
metabolic resource can be impaired or the metabolic 
resource is taken up by other organismic processes. 
While the reasons for a lowered dS/dt can differ, the 
effects on the mutation number probabilities are the 
same. 

Hypermutation under stress affects 
adaptation dynamics 

We already found that the minimally possible muta- 
tion rate increases when the metabolic resources uti- 
lized for the safekeeping of genetic information are 
limited. In our toy model for population adaption 
we can see the effects on adaptation dynamics by the 
comparison of two populations: Population (1): Both 
the growth rate G and the mutation rate (1 — p) are 
affected by the level of metabolic resources utilized. 
Population (2): Only the growth rate G is affected 
by the level of metabolic resources metabolized. 

Hypermutation effects in population com- 
parison: In Fig. [2] we see how stress-dependence 
of mutation rates influences adaptation dynamics in 
comparison to a conservative population with only 



growth rate depending on metabolic resource utiliza- 
tion. The first, most obvious effect is that on the 
speed of increase of the mean utilized metabolic re- 
sources dS/dt of both populations: In the example 
traces the stress-dependent mutation rate population 
adapts visibly faster than the conservative one with 
only growth rate differences. On a closer look, we can 
see the driving force a transient increase of mean mu- 
tation rate in the population with stress-dependent 
mutation rates. The pattern can be observed in many 
example traces, in the rough as well as mild type 
stress landscape scenario. 

C dependence of comparative adaptation 
dynamics: First, we compare for both stress land- 
scape types (mild and rough) the adaptation dynam- 
ics for different values of the parameter C, which 
characterizes the mutagenic influence from the en- 
vironment. In Fig. |3]we can see the effect of stress- 
dependent mutation rate increase. For the mild stress 
landscape, both populations (1) as well as (2) can be 
faster to adapt, dependent on the structure of the 
stress landscape. In the rough type stress landscape, 
faster adaptation of population (1) dominates, while 
for a some stress landscapes population (2) adapts 
faster than (1). For the rough type stress land- 
scape extreme acceleration effects are visible. For 
both stress landscape types the differences in adap- 
tation time course between both populations exhibit 
a strong dependence on parameter G . In the rough 
type stress landscape high C-values are associated 
with faster adaptation of population (1). In the mild 
type stress landscape population (2) mostly adapts 
faster at high G values. 

(C, D) dependence of comparative adapta- 
tion dynamics: Second, we compare for both stress 
landscape types (mild and rough) the adaptation 
dynamics for different values of the parameter du- 
plet (C,D). In Fig. U different typical regions show 
up: In the mild type stress landscape, the sam- 
ples describe a kind of "wedge" that expands from 
(C = 0.2, D = 0.2) and grows in D > direction 
for increasing G . Inside this wedge the (2) popu- 
lation exhibits faster adaptation dynamics. In the 
rest of the parameter space accelerated adaptation 
of the (1) population dominates. In the rough type 
stress landscape case no clear regions can be iden- 
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2: Example stress landscapes and resulting simulation traces: A) Example stress landscape of 
mild type m B) Example time course of adaptation for stress- induced hypermutation population (1) 
and conservative population (2) for mild type m stress landscape shown in A) C) Example stress 
landscape of rough type r D) Example time courses of adaptation for stress-induced hypermutation 
population (1) and conservative population (2) for rough type r stress landscape shown in C). For 
each stress landscape type (mild and rough) a stress landscape (^)g was drawn around a mean of 
D = 0.5. Colors, line styles and symbols: The ^ landscapes are shown as black vertical bars. In the 
same graph the mean of these bars is drawn as the horizontal dashed blue line. In the adaptation 
time courses, the solid (square symbols added for clarity in monochrome print) and the dashed black 
line (diamond symbols added for clarity in monochrome print) in the right graphs represent the 
mean fitness S' for population (1) and (2), respectively. The solid and the dashed red line (both 
without symbols for clarity in monochrome print) represent the mean per base pair mutation rate 
for population (1) and (2), respectively. Parameters used: D = 0.5, C = 1, A = 5, Euler time step 
0.1, simulation executed for 10000 time steps. 
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tificd, but mostly adaptation acceleration of popu- 
lation (1) is more frequent and the acceleration ef- 
fect is stronger. Sample (C,D) points showing faster 
acceleration of population (2) exhibit a trend to de- 
crease in occurrence in direction of increased D. Es- 
pecially in the rough type stress landscape, examples 
of hyper-adaptability occur, where up to 4-fold rel- 
ative acceleration of population (1) arise from the 
stress-dependence of mutation rate. 

Discussion 

The logic of the presented work is quite short and 
simple: 

A higher rate of primary mutations as well as a 
lowered ability to employ metabolic resources in muta- 
tion suppression increase the minimum effective mu- 
tation rate. This predicts transient mutation rate in- 
creases as a response to stress, which can accelerate 
adaptation. 

The logic is well-described by these two sentences, 
and the thermodynamical and quantitative analysis 
are just a more formal way of deriving and consol- 
idating this logic. Still, while the above logic is 
common fare amongst many biologists and medical 
researchers, the very specialists concerned with the 
topic have not actually reached a final common sense. 
Much experimental evidence indicates that the above 
logic applies [3S [TOl [28l [23 [26l El |2] and its conse- 
quences are discussed pljIS]. However, we have not 
come across a derivation of above logic from funda- 
mental physical axioms. Our work might add fun- 
damental physical credibility to the aforementioned 
experimental evidence of stress-induced hypermuta- 
tion and the resulting possibility of acceleration of 
adaption by transient mutation rate increase. 

We developed a purely thermodynamical, non- 
mechanistic analysis, which can by definition point 
out specific mechanisms. Instead, it can delineate 
the physical limitations, which biological information 
storage systems can not exceed. Experimental re- 
searchers find a whole range of mutagenic genome 
repair mechanisms and their resultant genetic mu- 
tation rates [m |37l [l^. A concept like the one we 
presented herein could possibly help to order some of 



this abundance of experimental knowledge. 

Specific selective scenario vs. piiysical con- 
straint: Hypermutation under stress and acceler- 
ated adaptation by transient mutation rate increase 
(Adaptive Mutation) demand an evolutionary expla- 
nation. 

First of all, it seems perfectly credible to use fre- 
quently changing living conditions of simple single 
cell organisms as a selective environment that favors 
adaptive agility. However, the explanation we sug- 
gest in this work is a physical constraint all cells are 
subject to, which does not require a specific selection 
scenario. It does not seem unreasonable, that mecha- 
nisms of genome repair and associated mutation sup- 
pression have developed to get close to the thermo- 
dynamical optimum. Both explanations, a specific 
selection scenario as well as the described physical 
constraint, lead to similar evolutionary development. 
However, latter explanation by a physical constraint 
needs weaker assumptions and applies more generally 
than one based on a specific selection scenario. Let 
us spell out the conditions which my analysis is based 
on: (1) The thermodynamical laws must hold, (2) an 
organism utilizes metabolic resources to maintain the 
integrity of its stored genetic information and (3) un- 
der stress this capability is compromised. This should 
most likely apply to all organisms and consequently 
the pattern of stress-induced hypermutation should 
be observable in many if not all of them. Note that 
a physical constraint and a specific selective scenario 
are in no way mutually exclusive. 

Discussion of findings in relation to otiier 
responses to compromised genetic integrity: 
A physical constraint as a reason for stress-induced 
hypermutation would be effective in all organisms. 
However, especially higher multi-cellular organisms 
can be expected to have developed counter-strategies 
to keep free of mutation effects in their constituent 
cellular organisms. Apoptosis or cellular senescence 
should be mentioned here. 

Gene repair and proofreading pathways are studied 
mostly in the context of DNA replication, see [21 ] l37 j . 
Differences in DNA replication also have pronounced 
effects on adaptation and most likely these effects and 
those from stress-induced hypermutation both influ- 
ence adaptation processes. 
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3: Comparison of adaptation dynamics for different C values: A) Mild type m stress landscape 
simulation. B) Rough type r stress landscape simulation. For each simulation run, the normalized 
half-time difference ST1/2 has been plotted vs. the sampled C value. ST1/2 > indicates a faster 
adaptation of the stress-induced hypermutation population (1), ST1/2 < indicates a faster adap- 
tation of the conservative population (2). For each condition (mild and rough) 500 simulation runs 
were executed. Parameters used: D = 0.5, N = 5, Euler time step 0.1, 500 simulations executed for 
20000 time steps each, only simulation runs with valid result shown. 
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Fig. 4: Comparison of adaptation dynamics for (C, D) sample points: A) Mild type m stress land- 
scape simulation. B) Rough type r stress landscape simulation. For each stress landscape type (mild 
and rough) 1000 simulations have been execiied. All those that reached sufficient growth during 
the simulation and for which |(^Ti/2| > 0.01 have been plotted into the above diagrams. Black lines 
represent examples in which the stress-induced hypermutation population (1) adapted faster, i.e. 

> 0. Red lines with a red circle represent a faster adaptation of the conservative population 
(2) adapting faster, i.e. 5Ti^2 < 0. Parameters used: N — 5, Euler time step 0.1, 1000 simulations 
executed for 15000 time steps, only simulation runs with valid results are shown. In the big regions 
without any sample points the simulation could not long enough to be properly evaluated. 



Discussion of findings in relation to other 
studies: We are aware of only two other theoretical 
studies that investigate the effects of stress-induced 
hypermutation on adaptation. [T] shows that in sta- 
ble, steady state populations stress-dependence in- 
creases genetic loads, as well as the cost of the main- 
tenance of sex. [13] investigates the influence an 
antiproportional relation between mutation rate and 
ornamentation has on sexual selection behavior and 
its effectiveness to avoid mates with germ-line mu- 
tations. In contrast, our work explores the physical 
credibility of stress-induced hypermutation as well its 
role in the time course of adaptation, which should 
prove especially interesting in medical scenarios of 
rapid pathogen evolution. 

While not directly concerned with the topic of this 
work, an abundance of theoretical literature exists on 
the role of mutation rates in evolution. Generally, 
mutation rates are assumed to be constant across 
genotypes. Some works do allow for mutation rate 
differences between genotvpes |18lBT| . Most interest- 
ingly, it has already been recognized, that popula- 
tions can exhibit accelerated adaptation by a tran- 
sient increase of mutator alleles, which exhibit the 
genetic trait of higher mutation rates [43]. In con- 
trast, a transient increase in mutation rate induced 
by stress is not a genetic trait, but simply an organ- 
ism's incapability to sustain safe storage of genetic 
information under stress conditions. Except [TJ 113] 
and our own approach presented herein. We know 
of no other theoretical work to investigate the role of 
such transient, stress- induced increases of mutation 
rates in adaptation. 

Also, the question if mutation rates are a se- 
lectable, genetic feature has been discussed at large, 
e.g see [121 |33l |3l [35l |6]. [7] gives a comprehen- 
sive overview how living conditions of an organism 
shape its mutation rate as a genetic trait. The 
bigger picture of recent studies on mutation rates, 
however, challenges the whole notion of a typical, 
mutation rate for a certain genotype, and instead 
suggests condition-dependent, dynamically changing 
mutation rates [141 122| . Our work might contribute 
fundamental credibility and insight to this notion 
of dynamically changing, condition-dependent muta- 
tion rates. 



A further related topic is that of the evolution 
of evolvabilitv[T5l [32] . and especially the causal- 
ity between variation and selection is a matter of 
concern. [16] As stated earlier, stress-induced hyper- 
mutation could indeed be a case of evolved evolv- 
ability, but the physical constraint we present in this 
work seems to be of more general applicability. The 
introduction of |16j gives an overview of experimen- 
tal results, theoretical work and discussion related to 
the evolution of evolvability. With regard to causal- 
ity, the interdependence of stress level and variation 
as per appearance of new mutants is reminiscent of a 
closed regulatory loop: A higher stress level induces 
a higher mutation rate, which allows swift adaptive 
mutation, which in turn relieves stress, and finally 
lowers the mutation rate again. Note, that this regu- 
latory circle arises from a fundamental physical con- 
straint, not as the result of a specific selective sce- 
nario. 

Another general class of works is concerned with 
the application of statistical mechanics methods to 
the study of the evolution on the DNA and RNA 
sequence level, e.g. see [1 [311 [24l |38] . Much of this 
work is of methodological value, and rests on assump- 
tions, which are challenged and rectified by works of 
the type of those mentioned shortly before. 

Experimental avenues: (1) Experimental stud- 
ies of the role of mutation rates in adaptation face a 
big conceptual problem: Mutations take effect along- 
side with other selective forces, e.g. growth and 
death rate differences. Our work suggest, that cells 
kept in non-growth but viable conditions still amass 
mutations in their DNA. Further, the probability 
and frequency of mutations should depend on the 
level of additionally imposed stress. Thus, an ex- 
perimenter could keep a population of cells on non- 
growth medium, while applying differential stress lev- 
els associated with different mutation probabilities 
and frequencies. Migration back to normal medium 
allows assessment of only the inflicted mutations, sep- 
arate of other adaptation effects. (2) A higher rate 
of primary mutations as well as a lowered ability to 
employ metabolic resources in mutation suppression 
can be caused by various stress types, the induction 
of hypermutation should be a common outcome for 
all of them. Thus, in the aforementioned experiment. 
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the exact type of stress should be irrelevant to the ob- 
servation of hypermutation. Further, different stress 
types should actually jointly increase the hypermuta- 
tion response, which suggests experiments with com- 
bined stress exposure. 

Medical implications: On an evolutionary scale, 
mutations are the source of development and biodi- 
versity. On the scale of the human individual, how- 
ever, mutations are mostly associated with undesir- 
able, often even critical consequences. A possible 
cause of cancer are genetic mutations and the re- 
sults presented in our work could indicate that in- 
creased stress in general increases genetic mutation 
rates and thereby the risk of cancer. Cancer resis- 
tances to chemotherapy as are supposedly driven by 
mutations. Our results indicate that chemotherapy 
could in fact induce hypermutation in cancer cells and 
thereby amplify their mutational adaptation. Simi- 
larly, our results support the view that antimicrobials 
do not only attack a pathogen but at sublethal doses 
also increase its adaptability by the stress inflicted 
upon the pathogen. |14j 

All mentioned phenomena are characterized by sur- 
prisingly swift adaptability of the respective cells, 
and the results from this work point towards stress- 
induced hypermutation as a possible common reason 
for this rapid adaptation. Firstly, this can deepen 
our understanding of the underlying causality, and 
secondly indicates, that therapeutic approaches used 
to disrupt the adaptability in one scenario, could pos- 
sibly be adapted to the other scenarios. 

Supplementary Materials 

MatLab scripts to produce the displayed simulation 
results. 

Acknowledgments 

The first idea was triggered by a discussion with Shi- 
Ian Mistry. Discussions with the following were es- 
sential: Siddharth Arora, Baboo Sabyasachi, David 
Hasselbach, Michael C. Mackey, Moises Santillan, 
Bartek Borek, Gregor Fussmann, Luke McNally, Pe- 
ter Richter and Isabel Gordo. All anonymous review- 



ers' efforts are deeply appreciated. 
Financial support from NSERC (Natural Sciences 
and Engineering Research Council of Canada), MI- 
TACS (Mathematics of Information Technology and 
Complex Systems) and a student and a traveling fel- 
lowship of the Studienstiftung des Deutschen Volkes. 
No funding body influenced the conception, execu- 
tion or writing of this work. 

References 

[1] AF Agrawal. Genetic loads under fitness- 
dependent mutation rates. J Evol Biol, 
15:10041010, 2002. 

[2] AF Agrawal and AD Wang. Increased trans- 
mission of mutations by low-condition females: 
evidence for condition-dependent DNA repair. 
PLoS Biol, 6(2):e30, 2008. 

[3] JB Andre and B Godelle. The evolution of muta- 
tion rate in finite asexual populations. Genetics, 
172:611626, 2006. 

[4] P Ao. Emerging of stochastic dynamical 
equalities and steady state thermodynamics 
from darwinian dynamics. Comm Theor Phys, 
49(5):1073-1090, 2008. 

[5] V Avila, D Chavarrias, E Sanchez, A Manrique, 
C Lopez-Fanjul, and A Garcia-Dorado. Increase 
of the spontaneous mutation rate in a long-term 
experiment with drosophila melanogaster. Ge- 
netics, 173(l):267-277, 2006. 

[6] CF Baer. Does mutation rate depend on itself? 
PLoS Biol, 6:e52, 2008. 

[7] CF Baer, MM Miyamoto, and DR Denver. Mu- 
tation rate variation in multicellular eukaryotes: 
causes and consequences. Nature Reviews Ge- 
netics, 8:619-631, 2007. 

[8] G Bell. Selection: The mechanism of evolution. 
Oxford University Press, 2 edition, 2008. 

[9] CK Bicbricher. Mutation, Gompetition and Se- 
lection as measured with Small RNA Molecules, 
chapter 4, pages 65-85. Academic Press, 1999. 



13 



[10] I Bjedov, O Tenaillon, B Gerard, V Souza, E De- 
namur, M Radman, F Taddei, and I Matic. 
Stressinduced mutagenesis in bacteria. Science, 
300:1404-1409, 2003. 

[11] L Brillouin. Science and Information Theory. 
Academic Press, 2 edition, 1962. 

[12] J Cairns, J Overbaugh, and S Miller. The origin 
of mutants. Nature, 335:142-45, 1988. 

[13] S Cotton. Condition-dependent mutation rates 
and sexual selection. Journal of Evolutionary 
Biology, 22:899-906, 2008. 

[14] A Couce and J Blazquez. Side effects of an- 
tibiotics on genetic variability. FEMS Microbiol 
Rev, 33:531-538, 2009. 

[15] A Crombach and P Hogeweg. Evolution of evolv- 
ability in gene regulatory netvirorks. PLoS Corn- 
put Biol, 4(7):el000112, 2008. 

[16] DJ Earl and MW Deem. Evolvability is a 
selectable trait. PNAS, 101(32):11531-11536, 
2004. 

[17] M Eigen and P Schuster. The hypercycle a prin- 
ciple of natural self-organization, part a emer- 
gence of the hypercycle. Naturwissenschaften, 
64:541-565, 1977. 

[18] Manfred Eigen. Selforganization of matter 
and the evolution of biological macromolecules. 
Naturwissenschaften, 10:465-523, 1971. 

[19] E Eisenberg, TL Hih, and Y-D Chen. Cross- 
bridge model of muscle contraction. Biophysical 
Journal, 29:195-227, 1980. 

[20] PL Foster. Stressinduced mutagenesis in bacte- 
ria. Crit Rev Biochem Mol Biol, 42(5):373-397, 
2007. 

[21] EC Friedberg. A brief history of the DNA repair 
field. Cell Res, 18:3-7, December 2008. 

[22] RS Galhardo, PJ Hastings, and SM Rosenberg. 
Mutation as a stress response and the regulation 
of evolvability. Critical Reviews in Biochemistry 
and Molecular Biology, 42:399-435, 2007. 



[23] RS Galhardo, PJ Hastings, and SM Rosenberg. 
Mutation as a stress response and the regula- 
tion of evolvability. Crit Rev Biochem Mol Biol, 
42(5):399-435, 2007. 

[24] T Garske and U Grimm. Maximum princi- 
ple and mutation thresholds for four-letter se- 
quence evolution. Journal of Statistical Mechan- 
ics: Theory and Experiment, page P07007, 2004. 

[25] H Ge and H Qian. The physical origins of 
entropy production, free energy dissipation and 
their mathematical representations. ArXiv, 
larXiv: 0911 .3984^2 [cond-mat/stat-mech] , 
November 2009. 

[26] S Goho and G Bell. Mild environmental stress 
elicits mutations affecting fitness in Chlamy- 
domonas. P Roy Soc London B, 267:123-129, 
2000. 

[27] BG Hall. Selection-induced mutations occur in 
yeast. Proc Natl Acad Sci, 89(10):4300-4303, 
1992. 

[28] PJ Hastings, A Slack, JF Petrosino, and 
SM Rosenberg. Adaptive amplification and 
point mutation are independent mechanisms: 
Evidence for various stress-inducible mutation 
mechanisms. PLoS Biol, 2(12):e399, 2004. 

[29] MN Hersh, RG Ponder, PJ Hastings, and 
SM Rosenberg. Adaptive muation and amplifica- 
tion in Escherichia coli: two pathways of genome 
adaptation under stress. Res Microbiol, 155:352- 
359, 2004. 

[30] TL Hill. Chemical Cycle Kinetics and Free En- 
ergy Transduction. Dover Publications, 2004. 

[31] GA Hoelzer, E Smith, and JW Pepper. On 
the logical relationship between natural selec- 
tion and self-organization. J Evol Biol, page 
17851794, 2006. 

[32] RE Lenski, JE Barrick, and C Ofria. Bal- 
ancing robustness and evolvability. PLoS Biol, 
4(12):e428, 2006. 



14 



[33] JI Lucas-Llcdo and M Lynch. Evolution of mu- 
tation rates: Phylogenomic analysis of the pho- 
tolyase/cryptochrome family. Mol Biol Evol, 
26(5):11431153, 2009. 

[34] SE Luria and M Dclbriick. Mutations of bac- 
teria from virus sensitivity to virus resistance. 
Genetics, 28(6):491-511, 1943. 

[35] AM Marcobal, DA Sela, YI Wolf, KS Makarova, 
and DA Mills. Role of hypermutability in the 
evolution of the genus oenococcus. J Bacteriol, 
190(2):564570, 2008. 

[36] A Meyerhans and J-P Vartanian. The Fidelity of 
Cellular and Viral Polymerases and its manip- 
ulation for Hypermutagenesis., chapter 5, pages 
87-114. Academic Press, 1999. 

[37] B Michel. After 30 years of study, the bacte- 
rial SOS response still surprises us. PLoS Biol, 
3:e255, July 2005. 

[38] J-M Park and MW Deem. Schwinger boson for- 
mulation and solution of the Crow-Kimura and 
Eigen models of quasispecies theory. Journal of 
Statistical Physics, 125(4) :975-1015, 2006. 

[39] L Perfeitol, L Fernandes, C Motal, and I Gordo. 
Adaptive mutations in bacteria: High rate and 
smah effects. Science, 317(5839):813-815, 2007. 

[40] I Prigoginc, G Nicolis, and A Babloyantz. 
Thermodynamics of evolution. Physics Today, 
25(12):38-44, 1972. 

[41] P Schuster and PF Stadler. Nature and Evolu- 
tion of Early Replicons, chapter 1, pages 1-24. 
Academic Press, 1999. 

[42] PD Sniegowski, PJ Gerrish, T Johnson, and 
A Shaver. The evolution of mutation rates: sep- 
arating causes from consequences. BioEssays, 
22:1057-1066, 2000. 

[43] F Taddei, M Radman, J Maynard-Smith, 
B Toupance, PH Gouyon, and B Godelle. The 
role of mutator alleles in adaptive evolution, na- 
ture, 387:700-702, 1997. 



[44] JS Wickcn. A thermodynamic theory of evolu- 
tion. J Theor Biol, 87:9-23, 1980. 

[45] JS Wicken. Evolution, Thermodynamics, and 
Information: Extending the Darwinian Pro- 
gram. Oxford Universtiy Press, 1987. 



15 



